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We present results of our quenched study of the B meson decay constant obtained with a parallel set of 
simulations with the Wilson and Clover actions at /3=5.9, 6.1 and 6.3. Systematic errors associated with the large 
6-quark mass are analyzed within the Fermilab non-relativistic formalism. As our best estimate in the continuum 
limit we obtain /s = 163±16 MeV and /s 3 =175±18 MeV with the Clover action. 
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1. Introduction 

A reliable determination of the B meson decay 
constant is a subject yet to be completed in lattice 
QCD. We have been pursuing this goal employing 
the relativistic formalism for heavy quark. Our 
results for the Wilson action has been reported 
in Refs. [§],§]. Since Lattice'96 we have carried 
out a parallel set of simulations with the O(a)- 
improved Clover action. We have analyzed the 
results for the two actions within the Fermilab 
non-relativistic formalism || with the view to un- 
derstand the systematic error due to a large value 
of 6-quark mass. In this article we report a sum- 
mary of results from the simulations and analyses. 

2. Simulation 

The parameters of our simulations are listed 
in Table [l]. For the clover coefficient we use 
the tadpole- modified one- loop value |j| given by 
Csw = p- 3 / 4 [l+0.199a y (l/a)] with P the average 
plaquette. Heavy quarks are simulated for 7 val- 
ues of the hopping parameter k covering the c and 



Table 1 

Simulation parameters. The lattice scale quoted 
is fixed by the string tension ^=427 MeV|Q. 



action 





5.9 


6.1 


6.3 




size 


16 a x40 


24 a x64 


32 a x80 




l/a (GeV) 


1.60(1) 


2.29(1) 


3.05(2) 




L (fm) 


2.0 


2.1 


2.1 


Wilson 


N con f 


150 


100 


100 


Clover 


N con f 


540 


200 


166 




Csw 


1.580 


1.525 


1.484 



* presented by S. Hashimoto 



b quark masses, and 4 values of k are employed 
for light quark. 

The heavy-light decay constant fp is extracted 
from the correlators of the axial vector current 
A4 and a smeared pseudoscalar density P s (x) = 
^2p4>{\r\)Q(x + r)~fsq(x) on the Coulomb gauge 
fixed gluon configurations, where </>(|r|) is the 
pseudoscalar wave function measured for each 
heavy and light quark masses. 

A new perturbative ingredient in our work is 
the recent one- loop result Q for the pole (mf) 
and kinetic (m^) masses of heavy quark and the 
renormalization factor 2^(omQ) of the axial vec- 
tor current for finite bare heavy quark mass tuq 
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Figure 1. $p/ v / ct 3 ^ 2 as a function of *Ja/mp. 
Filled (open) symbols are Clover (Wilson) results. 
Circles, squares and diamonds correspond to val- 
ues at /3=5.9, 6.1 and 6.3. 

(see Refs. for previous calculation for the 

Wilson case). Effects of finite values of airiQ in 
Za is significant, reducing fp by 5-2% for the 
Wilson action and increasing it by a similar mag- 
nitude for the Clover case compared to the value 
obtained with Zji(ainQ — 0). 

We define the heavy-light meson mass by|]J^] 
TTi p = rripi+m® — with mp\ the measured 
meson pole mass and the one-loop perturbative 
result |j applied for m^ — m^. This definition 
does not have the problem of the measured kinetic 
mass that the b quark mass can not be deter- 
mined consistently from heavy-light and heavy- 
heavy mesons |lQ,nfi,p[ . 



3. Results 

We plot the quantity defined by $(rap) = 
(a s (mp)/a s (mp)) 2 ^ fpy/mp in Fig. [I] as a func- 
tion of 1/mp for the the Clover(filled symbols) 
and Wilson (open symbols) actions. The light 
quark mass is linearly extrapolated to the chi- 
ral limit, and a s ()i) is calculated with the stan- 
dard 2-loop definition with Aqcd = 295 MeV. 
We normalize the results by the string tension 
ctQ since we primarily wish to examine the ques- 
tion of large-amQ errors in this figure. Vertical 
lines indicate the position of B and D mesons for 
%fo = 427 MeV. Data points at 1/mp = are 
the static results p2[, to which our results for the 
same set of (3 values converge. 



Figure 2. Continuum extrapolation of fp (cir- 
cles) and fp (squares). Filled (open) symbols are 
Clover (Wilson) results. 

In Fig. H we plot the continuum extrapolation 
of fp and fp. The Wilson results exhibit a scal- 
ing violation of 11-5% in our range of lattice spac- 
ing a -1 1.6-3 GeV, while the Clover results 
show a significantly reduced variation of 4-2%. 
These magnitudes are common to fp and fp. 
Furthermore the continuum values obtained with 
the two actions by a linear extrapolation agree 
within the statistical error of about 5%. 

We emphasize that this agreement does not 
necessarily mean that systematic errors due to 
heavy quark are negligibly small. In the non- 
relativistic interpretation, the equivalent Hamil- 
tonian for Wilson-type actions has the form 



H = Q 



mi 



D 2 

2my 



ia B 
2m B 



0(l/m 2 Q ) 



Q- (1) 



For the Wilson action for which mp ^ rri2, the 
leading error in the decay constant due to heavy 
quark is 0((cp-l)kQCp/m Q ) with cp = m 2 /mp. 
For the B meson an examination cp at the tree 
level shows that a linear extrapolation of cp from 
our range of a^ 1 leads to a value |cp — 1| w0.4 at 
a=0. We should therefore allow an 0(3%) error 
unremoved in the continuum limit where we used 
Aqcd=0.3 GeV. The same magnitude of error 
also remains for fp. 

There are two more sources of systematic error 
we need to consider. One is mQ-independent scal- 
ing violation of O(clAqcp), which we estimate to 
be 10% at our smallest a -1 «3 GeV. The other 
is 0(a 2 s ) uncertainty due to the use of one-loop 
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Figure 3. Ratio of lattice scale obtained from 
m p (circles) and from /„. (squares) to that from 
the string tension. Filled (open) symbols are for 
Clover (Wilson) action. 

value for Za, which is 0(4%) with a s « 0.2. 
Adding all the errors by quadrature leads to a 
combined systematic error of 0(11%) in the de- 
cay constant for the Wilson case. 

For the Clover action for which ms = vni 
to 0(a s ), the large-omQ errors have the form 
0(a s AQCD/m Q ) and 0(K 2 QCD /ni 2 Q ). We esti- 
mate their magnitude to be 0(1%) incorporating 
the effect of coefficients, similar to Cb, that van- 
ish in the continuum limit. The scaling violation 
errors are 0{a s ah.Q C D, a2 ^ 2 qcd) which are small 
at 0(2%). With the 2-loop error of 0(4%) from 
Za and an additional error of 0(2%) arising from 
the field rotation ignored in the present calcula- 
tion, the combined systematic error amounts to 
0(5%) for the Clover case. 

So far we have used the string tension a to set 
the scale. In Fig. || we plot the ratio of a -1 ob- 
tained with m p and f n to that with a. We use 
the variation of the ratio to estimate the uncer- 
tainty in setting the scale, which we take to be 
10% for the Wilson case and 5% for the Clover 
case. This uncertainty includes possible quench- 
ing error as the ratio need not converge to unity 
in the continuum limit. 

Our final result for the decay constant is tab- 
ulated in Table |[ To obtain the values we take 
the continuum extrapolation of fp/y/a and con- 
vert it fp/m,p with the value of \fa/m p at a=0 
in Fig. |3|. A direct continuum extrapolation of 
fp/mp yields consistent results. The quoted er- 
rors are statistical, systematic and due to scale 



Table 2 

Results for the decay constant in MeV unit. 





Wilson 


Clover 


fB 


140(11) (15) (24) 


163(9)(8)(11) 


f Bs 


159(10)(17)(27) 


175(9)(9)(13) 


Id 


163(13)(18)(28) 


184(9)(9)(12) 


fD. 


180(11) (20) (31) 


203(9) (10) (14) 



setting as estimated above in this order. We 
take the result for the Clover action to be our 
best estimate primarily because scaling violation 
is smaller and also since the statistical ensemble 
is larger compared to the Wilson action. Com- 
bining errors by quadrature we obtain the results 
quoted in the abstract. 
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